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SUMMARY 

Computer curve fitting is used in mul ti wave length pyrometry to measure 
the temperature of a surface in the presence of reflected radiation by decom- 
posing its radiation spectrum. Computer-simulated spectra (at a surface tem- 
perature of 1000 K; in the wavelength region 0.3 to 20 pm; with a reflected 
radiation-source temperature of 700 to 2500 K and reflector emissivity from 
0.1 to 0.9) were generated and decomposed. This method of pyrometry determined 
the surface temperatures under these conditions to within 5 percent. The prac- 
ticability of the method was further demonstrated by the successful analysis ot 
a related problem - decomposition of the real spectrum of an infrared source 
containing two emitters to determine their temperatures. 


INTRODUCTION 

In pyrometry, the temperature of a surface is determined by measuring the 
emitted radiation. The emitted radiation obeys Planck's blackbody distribution 
law. In conventional pyrometers (the so-called single-color pyrometers), the 
radiation is measured in a narrow wavelength band, usually in the visible or 
near-infrared bands (in the region of 0.5 to 1.0 pm). This method of temper- 
ature measurement requires some knowledge of the surface emissivity. Two-color 
pyrometers measure the emitted radiation in two wavelength bands, and by assum 
ing that the ratio of emissivities in the two bands is known, the surface tem- 
perature is determined without knowing its emissivity. In general, as long as 
the emissivity does not change rapidly with wavelength, accurate temperature 
measurements are obtained. However, both methods are susceptible to errors 
caused by reflected radiation. Reflected radiation originates from sources 
other than the surface under consideration and reaches the detector through 
reflection from the surface. Ceramics are materials from which future advanced 
propulsion systems will be built. This class of material typically has low 
emissivity (0.15 to 0.8); hence, there is high reflectivity in the wavelength 
band where conventional pyrometers operate (ref. 1). Consequently, a potential 
for large error exists when pyrometry is employed to measure the temperature ot 
ceramics. To correct for the reflected radiation, we propose the use of multi- 
wavelength pyrometry that measures temperature by curve fitting the measured 
radiation spectrum. This method has been used to simultaneously measure tem- 
perature and spectral emissivity with accuracy up to 1 percent (refs. 2 to 4) 
and to measure temperature distribution with a known spectral emissivity 
(ref. 5). 


REFLECTED RADIATION 

Figure 1 shows a pyrometry arrangement in which both the emitted and 
reflected radiations reach the detector. One can see that, in addition to the 


radiation (signal) originating from the measured surface at temperature T e , 
there is an undesirable component originating from an extraneous source at tem- 
perature T r . The arrangement shows that the radiant flux detected by a pyrom- 
eter can be represented mathematically as 

<t>X = 4>x,e + 4>X,r (1) 

where <j>x is the detected flux decomposed into two components: <(>x e , caused 

by emission from the surface under consideration and 4>x >r , caused 6y reflec- 
tion from the extraneous source. As a result, depending ’on the magnitude of 
this flux, the final temperature measurement contains unavoidable errors. 
Neither the one-color nor the two-color pyrometry is capable of accounting for 
and, hence, eliminating this source of error. 

By application of Kirchoff's law, and the Planck distribution of blackbody 
radiation Lx(T), the two components are expressed as 

, e “ ae e(^)L\(Te) 

^X , r = a bo e (X)e r (X)Lx(T r ) 

where 

c l l 

f f t ) = _t I 

X* ^5 exp(c 2 /XT) - 1 

In equations (2) and (3), the quantity a is a geometric constant of the py- 
rometer, and b is a geometric constant of the experimental setup containing 
quantities such as the view factors (solid angle). The radiation constants 
ci and C 2 are given by 

c t = 2c 2 h (4) 

c 2 = ch/k (5) 

where c = 3xl0 8 m/sec is the velocity of light, h - 6.6xl0~ 34 J-sec is the 
Planck constant , k = 1.38xl0~ 23 J/K is the Boltzmann constant, and e e (X), is 
the emissivity of the surface, which is generally a function of wavelength and 
may also be temperature dependent. The quantity e r (X) is the emissivity of 
the extraneous source, and for a blackbody, e r (X) = 1. The quantity a e (X) is 
the reflectance of the target surface and is equal to (1 - e e (X)), if it has no 
specular component. Multiple reflection between the two surfaces is neglected. 
Curves of the Planck radiation distribution for a blackbody are shown in fig- 
ure 2 as a function of wavelength for different temperatures. A temperature T 
can be calculated from the relative spectral distribution near the peak because 
the wavelength at the peak varies inversely with T (Wien's displacement law). 

It is thus possible to find the two temperatures T e and T r in figure 1 
because of their different spectral distributions. 


PYROMETRY METHODS 

One- and Two-Color Pyrometry and Polaradiometry 

Temperature measurement by pyrometry is based on detecting and analyzing 
the radiation of equation (1). In both one- and two-color pyrometry, the 


( 2 ) 

(3) 
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reflected radiation is ignored or considered small compared to the emitted 
radiation. This is not a bad approximation for two-color pyrome try, if the 
extraneous radiation sources have temperatures not much different from the 
target surface. 

In the single-color pyrometer, the radiation emitted in a narrow wave- 
length band AX centered at \i is detected by a suitable detector system. 
The detected flux <t> is equal to the integral obtained by integrating equa- 
tion (2) over this wavelength region: 


H\) = 


p X 1 +AX/2 


e e(X) L x( T e )dX 


J X 1 -AX/2 


( 6 ) 


When AX is small, it can be assumed that the integrand in equation (6) is 
constant and that the integral is approximated by AX. Then the detector 
signal is 


S 1 = k i aAXe e (V x 5 exp(c 2 /X 1 T e ) - 1 (7) 

where ki is the detector sensitivity factor. By knowing the emissivity of the 
body, we can determine its temperature from the detector signal. 

If a second narrow wavelength band AX centered around X 2 is selected, 
the detected signal S 2 is 


s 2 “ k 2 aAXe e (X 2 ) x 5 exp(c 2 /X 2 T e ) - 1 (8) 

The ratio of equations (7) and (8) contains the ratio of the emissivi t ies , 
which is unity if they are the same. This is the two-color pyrometry method. 
Neither method addresses the issue of reflected radiation. If it is not negli- 
gible, the resulting error can be substantial. 

Murray (ref. 6) demonstrated an ingenious method for obtaining the temper- 
ature of a surface by explicitly employing the reflected radiation from a known 
external source. By carefully analyzing the emitted and reflected radiations, 
explicitly taking into consideration their polarizations (polarization paral lei 
and perpendicular to the incident plane), and applying Kirchoff's law in a 
narrow band around A, Murray obtained an expression for the signal I caused 
by the combined radiation. 



In his expression, ep and e s are the sample emittances for radiation 
polarized parallel and perpendicular to the plane of incidence; Eg is the 
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square of the amplitude of the radiation of wavelength X from the blackbody 
with temperature T; the subscripts R and S refer to the reflecting surface 
and the reference source; and 9 is the polarizer analyzing angle. 

This polaradiometer signal contained a constant dc component and a sinu- 
soidal ac component proportional to the difference between emissivi ties and 
the difference between the emitted and reflected radiations. By substituting 


e e (X) ‘ < e p * e s > 

(10) 

( 1/4 )E bs - exp^/xj^j _ ! 

(U) 

( 1/4 ) e br » exp^/xTj _ i 

(12) 


into the dc component, the quantities in the dc component are identical to the 
results represented by equations (1) to (3). 


Pyrometry by Curve Fitting 


Unlike Murray's polar-pyrometry which depended on knowledge of the exter- 
nal blackbody temperature to measure the temperature of the surface, the pro- 
posed pyrometry technique will be able to correct for the interfering effect of 
reflected radiation and to measure the temperatures of the surface and the 
extraneous source simultaneously. In polar-pyrometry, as in the one-color 
pyrometry, radiation is detected in only one, narrow, wavelength band. The 
proposed technique calls for measuring the intensity of the radiation in the 
wavelength region on both sides of the peak in the blackbody curve not - just 
at a narrow band, but at many wavelengths. 

A spectral radiometer is well suited to obtaining a spectrum represented 
by equation (1). It is a spectrometer which accurately measures the radiation 
emitted by a body as a function of wavelength. It has high spectral-resolution 
capability, with spectral steps as fine as 5 A. Once a spectrum is obtained, 
statistical techniques in spectroscopy and other fields of research can be used 
to extract the parameters such as a, b, T e , and T r by fitting the mathemati- 
cal expression of equation (1) (with eqs. (2) and (3) substituted into it) to 
the experimental radiation spectrum. In curve fitting, an initial set of 
approximate values for the unknown parameters is assumed. From these initial 
values, the expression given by equation (1) is evaluated. The residuals, 
which are the differences between the calculated and the experimentally 
observed quantities, are obtained. Next, the values of the unknown parameters 
are varied according to some mathematical algorithm, and a new set of residuals 
obtained. This is performed for many iterations so that the sum of the squares 
of the residuals gets progressively smaller, and the set of parameter values 
which produces the least squares will be the best estimate. In this way, the 
experimentally obtained spectrum has, in fact, been decomposed into its 
emitting and reflecting terms (i.e., the temperatures are simultaneously 
determined) . 
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COMPUTER EXPERIMENTATION 


Simulated composite spectra were generated mathematically according to 
equations (1) to (3). These spectra consisted of radiation from a surface 
with emissivity e e (X) = constant and temperature T e , and reflected radiation 
from a blackbody with emissivity e r (X) = 1 and temperature T r . The chosen 
temperature T e was 1000 K and T r values were 700, 1500, 2000, and 2500 K. 
Only three significant figures were retained by this calculation. This numeri- 
cal roundoff is equivalent to the introduction of some error into the data. 

The wavelength region of the simulated spectrum is from 0.3 to 20 pm repre- 
sented by 37 equally spaced channels. Values of e e vary from 0.1 to 0.9. 

The temperatures of the two surfaces and the emissivity e e were deter- 
mined from the spectra. These were decomposed on a personal computer by 
running a commercial, scientific, statistics package (with least-squares, 
curve-fitting capabilities) called RSI. Other similar computer software could 
also have been used. 

Curve fitting decomposes a spectrum into its constituent Planck functions 
by simultaneous extraction of the two temperatures. The computer program also 
yields the standard errors of these quantities resulting from the statistics of 
fitting the data to the functional relationship assumed by equations (1) to 
(3). The reported significance level is 0.0001. In the present computer 
experimentation, we define the error as the deviation of the fitted temperature 
from 1000 K. These errors are of the same order-of-magni tude as those esti- 
mated statistically by the computer program. The variation of these errors 
with temperature and emissivity is shown in figure 3. It shows that as the 
emissivity decreases, the fractional error of the result increases rapidly. 

The results indicate that for emissivity larger than 0.3, there is sufficient 
accuracy in the temperature so determined. 

By increasing the number of channels in the wavelength region represented 
by the spectrum, the fractional error can be reduced. The case of T e - 1000 If, 
T r - 1500 K, e e = 0.1, and bandwidth 0.3 to 20 um is shown in figure 4. When 
the number of channels is increased from 37 to 193, 394, 462 and 924 (corre- 
sponding to channel spacing of 0.5, 0.2, 0.1, 0.05, and 0.02 pm), the frac- 
tional error decreases from 0.05 to less than 0.01. However, this increase in 
accuracy is accompanied by a corresponding increase in computation time. 

The fractional error increases when the spectral bandwidth describing the 
spectrum is reduced. In addition to the 0.3- to 20-pm wavelength region, the 
following wavelength regions were also investigated (fig. 2): (1) 1.8 to 6 pm 

(corresponding to detecting radiation greater than 30 percent of the 1000 K 
peak), (2) 1.2 to 10 pm (corresponding to detecting radiation greater than 
10 percent of the 1000 K radiation peak), and (3) 1.0 to 15 pm (corresponding 
to detecting radiation greater than 2 percent of the 1000 K peak). The results 
in figure 4 show that the bandwidth should be at least 1.2 to 10 pm (a wave- 
length ratio of 8) . 

A channel spacing of 0.05 pm seems to be an optimal value. The fractional 
error for various extraneous source temperatures, as a function of emissivity 
with this channel spacing, is shown in figure 5. It can be seen that except 
for an extraneous source temperature of 2500 K and surface emissivity 0.1, the 
fractional error for all the other situations is less than 0.04. The results 
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indicate that for this optimized channel spacing of 0.05 pm, there is suffi- 
cient accuracy in the temperature determined when the emissivity is larger than 
0.2, compared to the necessary value of emissivity larger than 0.3 mentioned 
previously when the channel spacing was 0.5 pm. 

It is informative to compare the flux $ e from the surface to the flux 
<J> r from the reflected extraneous source. Integration of equations (2) and (3) 
over the range of 0.3 to 20 pm is approximated by integration from zero to 
infinity. The resultant flux follows the T 4 Stefan-Bol tzmann law. The flux 
ratio 4> e /<t> r is given by 


e e^> 

ba e (X) 
e r l ’ 



(13) 


For the case e ? (X) = 0.1, e r - 1, o e - 1 - e e (X), T e = 1000 K, T r = 1500 K, 
and b - 1 (b is given by the ratio of the solid angle subtended by the 1500 K 
body and 2ir), this flux ratio is 0.022. Thus the temperature can be deter- 
mined in the presence of high levels of reflecting flux. 


INSTRUMENTATION 

A spectral radiometer is used to obtain a complete spectrum of the body 
whose temperature we want to measure. Several commercial units exist that 
operate from about 0.3 pm to as far as 20 pm. These instruments are completely 
automated by microprocessor control. The resultant spectrum is in digitized 
numerical form and can be easily transferred to a computer system such as the 
IBM, IBM PC, and the VAX for processing and analysis. 

It may take about 30 sec to accumulate a spectrum, but for development 
purposes, this is quite adequate. The collected spectrum is then analyzed by 
a statistical program. However, in many applications, the 30-sec data- 
acquisition time may be too slow. In that situation, an optical multichannel 
analyzer, which accumulates the whole spectrum instantly, can be used. 

Presently, most optical multichannel analyzers use array silicon detec- 
tors, which are sensitive from about 0.2 to 1 pm. The spectra of transient 
events are routinely collected by workers in spectroscopy. For pyrometry pur- 
poses, in order to detect the whole spectrum, operation up to at least 2.5 pm 
is required. Fortunately, newer detectors, such as InGaAs and PtSi, are 
becoming available in array packages. Thus, the data acquisition can be 
instantly accomplished. 

The temperature extraction is accomplished by computer software. In the 
preliminary results documented here, the data transfer and interfacing between 
the computer programs are performed manually. In a fully developed system, a 
new-generat ion , high-speed, personal computer will be incorporated into the 
data acquisition with optimized software. Temperature output in a matter of 
seconds is no doubt possible. 
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SOME EXAMPLES 


The mul t i wave length pyrometry technique discussed in the "Pyrometry 
Methods" section was applied to analyze the real spectrum of two radiators. 

The problem is mathematically identical to the simulation described in the 
"Computer Experimentation” section. 

There are two, commercial, Oriel, infrared sources that have data sheets 
providing the spectral irradiance from 1 to 20 pm (Oriel Corporation Catalog, 
Vol. II, Light Sources, Monochromators, Detector Systems, 1989.): the 22-W, 

1550 K, 6575 ceramic element and the 100-W, 1000 K, 6363 infrared emitter. 

The data sheets provide enough information to determine temperatures of these 
sources (fig. 6). The data for these infrared sources were fitted to the sin- 
gle Planck curve (eq. (2)) and constant emissivity was assumed to obtain tem- 
peratures, of 976 K (within 2.5 percent error) for the 6363 ceramic element and 
1375 K for the 6575 infrared emitter (with an unacceptable error of almost 
12 percent). Figure 7 shows the fitted curves and data for the two infrared 
sources . 

The agreement for the 6363 source was very good; whereas, the agreement 
for the 6575 was certainly unacceptable. An examination of the graph in the 
data sheet revealed some structures in the spectrum between 3 and 6 pm; the 
cause was suspected to be the variability in the 6575 emissivity with wave- 
length. Indeed, the manufacturer has provided, in the form of a broadband 
curve, the emissivity of the 6575 ceramic element as a function of wavelength 
(fig. 8). The emissivity data in that graph was fitted empirically to a 
Gaussian relation 



" 2a 2 ) 


(14) 


producing a maximum at X 0 = 14 pm and a = 7.75 pm. 

This expression for the emissivity was then inserted explicitly into the 
Planck formula in equation (2), and the data was fitted accordingly. The tem- 
perature obtained is 1571 K, with an error of less than 1.5 percent from the 
nominal 1550 K. The fitted results, shown in figure 9, are an improvement from 
those shown in figure 7; nevertheless, in the region between 4 and 14 pm, the 
experimental data in figure 9 systematically exceed the calculated value 
obtained by curve fitting. 

To explain this, we assume the presence of another radiator peaking at 
about 5 pm buried under the original curve. When the radiation from the second 
radiator is treated as if it were reflected from a surface described in the 
section, "Computer Experimentation," the mathematical expression describing the 
radiation from these two radiators can be immediately written as 

L - L el + L e 2 ( 15 ) 


or 


L 


C 1 ! c 3 i 

e l^ ,5 exp(c„/XT) - 1 + e 2^ -,5 exp(c„/XT„) 


(16) 
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where T and T 2 are the temperatures of the two radiators. The constants 
ci and C 3 now absorb the first radiation constant and the surface area of 
the emitting sources. Curve fitting produced c^ = 1602, T - 777 K, C 3 = 70, 
and T 2 = 1642 K. The final fitting, shown in figure 10, is in excellent 
agreement with the data. 

The decomposed spectrum, shown in figure 11 , indicates that the infrared 
source 6575 actually contains two heating elements: one operating at a temper- 

ature of 1642 K and the other operating at 777 K. The emitting surface areas 
are approximately 20 to 1 . The manufacturer confirmed that this infrared 
source is constructed with platinum wire wrappings buried in a cylinder of 
refractory material. During operation, there are effectively two heating ele- 
ments exposed, and the ratio of their radii is about 10 to 1 . 


CONCLUSION 

Computer curve fitting has been used in multiwavelength pyrometry to meas- 
ure the temperature of a surface in the presence of reflected radiation. The 
technique was successfully applied to a computer-generated composite spectrum 
as well as real commercial sources. The ability of multiwavelength pyrometry 
to deconvolute the spectrum demonstrates the applicability of this technique to 
analysis of the complex spectrum of low-emissivi ty materials such as ceramics. 
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■ Pyrometer 
Figure 1 . - Typical pyrometry geometry. 



Figure 2. - Blackbody radiation curves. Wavelength 
Bands that included radiations greater than 2, 10, 
and 30 percent of the 1000 K peak are shown. 



Figure 3. - Fractional error in surface temperature in 
presence of a blackbody as a function of emissivity. 
Surface temperature, 1000 K; channel spacing, 

0.5 jam; and spectral bandwidth, 0.3 to 20 jim. 



Figure 4. - Fractional error of surface temperature in 
presence of a blackbody as a function of channel 
spacing in wavelength region. Surface temperature, 
1000 K; blackbody temperature, 1500 K; and emis- 
sivity of the surface, e e = 0. 1 . 
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Temperature, 



Figure 5, - Fractional error in surface temperature in 
presence of a blackbody as a function of emissivity. 
Surface temperature, 1000 K; spectrum channel 
spacing, 0.05 pm; and spectral bandwidth, 0.3 to 
20 pm. 



Figure 6. - Spectrum of two infrared sources. Typical 
spectral irradiance from bare element per 10-mm 2 
area. (Oriel Corporation Catalog, Vol. II, Light 
Sources, Monochromators, Detector System, 1989.) 




Figure 7. - Spectra of two infrared sources • fitted to 
Planck curves. 


Figure 8. - Spectral emissivity of 6575 ceramic-element 
infrared heater. (Oriel Corporation Catalog, Vol. II, 

Light Sources, Monochromators, Detector System, 1989.) 
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Figure 9. - Spectrum of 6575 ceramic- element lamp 
fitted to a Planck curve of temperature 1571 K and 
incorporating the spectral emissivity data. 



Figure 10. - Spectrum of 6575 ceramic-element lamp 
fitted to two Planck curves at temperatures of 1642 
and 777 K and incorporating the spectral emissivity 
data. 



Figure 11.- Deconvolution of the spectrum of the 
6575 ceramic-element infrared heater into the 
spectra of two Planck functions and incorporating 
experimental data. 
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